Sudden variations in the speed of sound during inflation: 
features in the power spectrum and bispectrum 
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We employ the formalism ol the effective field theory ol inflation to study the effects of a sudden 
change in the speed of sound of the inflationary perturbations. Such an event generates a feature 
with high frequency oscillations both in the two- and in the three-point functions of the curvature 
fluctuations. We study, at first order in the magnitude of the change of the speed of sound, the 
dependence of the power spectrum and of the bispectrum on the duration of the change. In the limit 
of a very short duration, the oscillations in the power spectrum persist up to very large momenta 
and the amplitude of the feature in the bispectrum diverges while its location moves to increasing 
momenta. 
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I. INTRODUCTION 

Inflation, in its simplest realizations, is a period of 
rather uneventful quasi-de Sitter expansion of the Uni- 
verse. As a consequence, it typically predicts a quasi- 
scale invariant, featureless spectrum of scalar and tensor 
metric perturbations. It is however possible that the in- 
flating Universe did experience some traumatic events 
that left a permanent mark on the perturbations we cur- 
rently observe in the Cosmic Microwave Background ra- 
diation. 

Various examples of such traumatic events have been 
discussed in the past. The simplest possibility is that 
they originate from a feature, possibly associated to a 
phase transition [1], in the shape of the inflationary po- 
tential [2-6] . Another possibility is that the slowly rolling 
inflaton induces an explosive production of quanta of 
some other degree of freedom [7-9] that affect the met- 
ric perturbations. In general these phenomena lead to 
high frequency oscillations in the primordial power spec- 
trum about the scales that left the horizon at the time 
of the event. Several groups - see for instance [10-12] - 
have found hints of the presence of such oscillations in 
the CMB data. 

In the present paper we study a different kind of sud- 
den event: a change in the speed of sound of the infla- 
ton perturbations. The speed of sound c s of the inflaton 
can be different from the speed of light in models with 
nonstandard kinetic terms [13], and in particular in DBI 
inflation [14]. It is therefore natural to ask what observ- 
able effects can occur when the speed of sound suddenly 
changed during inflation. In the past, this situation has 
been considered in [15] and in [16], which studied the 
effect of a jump in c s on the two-point function of the in- 
flaton. Reference [15] did also estimate the effect of such 
a jump on the bispectrum. 
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In our work we provide a detailed study of how a rapid 
but not instantaneous change in c s affects both the power 
spectrum and the bispectrum of the curvature pertur- 
bations. The Lagrangian provided by the effective field 
theory of the inflaton fluctuations [17] serves our purpose 
well, in that it allows us to determine the effect of the 
jump in the speed of sound in a model-independent way. 
Differently from what happens in most of the literature, 
we study the system of [17] without assuming that the 
background parameters evolve adiabatically. 

Since we work in the context of [17], we are not in- 
terested in the details of the full theory producing the 
sudden change in c s , and will not worry about the is- 
sue of its technical naturalness. Let us just note that 
reference [15] discusses an explicit example where brane 
inflation leads to a jump in c s . Reference [16] contains 
a couple of additional examples. And here we mention 
one more possibility: In general, the speed of sound in 
DBI inflation [14] depends on the expectation value of 
several fields [18]. The rolling inflaton can lead, a la [7], 
to the nonperturbative production of matter in the sec- 
tor determining the value of c s . If this occurs, then the 
speed of sound will momentarily change its value before 
the nonperturbatively produced particles dilute away. 

Let us now summarize our main results. The back- 
ground dynamics producing a change in the speed of 
sound will also induce, in general, a change in the slow 
roll parameters e and r\. In order to disentangle the dif- 
ferent effects, however, we will consider a scenario where 
e and r\ are kept constant and only c s changes. Now, as 
long as both the speed of sound and the slow roll param- 
eters are constant, the amplitude of the power spectrum 
is inversely proportional to c s and to e. The spectral in- 
dex, on the other hand, depends on e and rj, but not on 
c s . As a consequence, a sharp jump in c s will induce a 
jump in the amplitude of the power spectrum without 
affecting its tilt. Our analytical and numerical work will 
confirm this expectation. This property might be used 
to discriminate models with a sudden change in the in- 
flationary potential from those with a sudden change in 
c s . One might even be tempted to use a sudden change 
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in the speed of sound to explain the small amplitude of 
the observed CMB fluctuations at the largest multipolcs. 

Our system depends on two parameters - the change 
Ac s in the speed of sound and the duration At of such 
a change. As long as At is not vanishing, the two-point 
function shows an oscillating pattern that extends over a 
finite range of momenta. The duration of the oscillations 
becomes infinite in the limit At — > 0, in agreement with 
the results of [15, 16], reflecting a UV artifact that orig- 
inates from the discontinuous behavior of c s . The three- 
point function or /nl, when evaluated in the equilateral 
configuration, also oscillates, with a maximum amplitude 
scaling as Ac s Tq/At 2 at momentum ~ At -1 , where To 
is the (conformal) time when the change in c s occurs 
(the amplitude of the bispectrum in the squeezed limit 
scales as Ac s |to|/At, and agrees with the consistency re- 
lation [19]). Therefore a jump that is rapid enough may 
cause large nongaussianity at large momenta. The be- 
havior of the bispectrum turns out to be similar to that 
produced by a sudden change in the inflaton potential re- 
cently studied in [6] , and future surveys might be able to 
detect it if the feature appears at sufficiently large scales. 

The plan of the paper is the following. In section II 
we apply the gauge invariant formulation of cosmological 
perturbations to the setup of [17], and derive a quadratic 
action for the gauge invariant variable v describing the 
scalar perturbation. Then we model the jump in c s and 
solve the equation of motion for v, to obtain the expres- 
sion of the power spectrum. In section III we take the de- 
coupling limit, study its regime of validity and calculate, 
in this limit, the equilateral and the squeezed bispectra 
for our system. Finally, in section IV we summarize our 
results and discuss future prospects. 



the action (1) is straightforward. At first order in per- 
turbations we obtain the following background equations 



c + A-3M£H z =0, 

c - A + M 2 (3iJ 2 + 2H) = , 



(3) 



with an overdot implying a differentiation with respect 
to the coordinate time t and H = a/a. 

Before we get to the second order, we write <j> and B 
in terms of the other field variables by perturbatively 
solving the lapse and shift constraints SS/SN = and 
SS/SN* = 0: 



eHa 



H 



(4) 



A£> 

a z a 2 H 



where we have defined 



H 



2 Mi 



e ~ h 2 ' mi h 2 ' 

and where the speed of sound is given by 

e 



e + fi 



(5) 



(6) 



With (4) plugged in, the quadratic action depends only 
on ip, E and a. The terms in E, however, can be dropped 
because they are total derivatives. Finally, we introduce 
the gauge invariant variable 



v = z (ip + Ha) , 



(7) 



II. GAUGE INVARIANT FORMULATION AT 
QUADRATIC ORDER 

A. Action 

We employ the effective field theory of inflation pro- 
posed by [17]. Since we are interested only in changes in 
the speed of sound, we neglect the terms in M 3 , M 2 and 
M 3 , etc. of [17] and start with the following action 



with 



M 2 

-A(t + a) + ^M 2 (t + <j)\l+g m ) 2 



(1) 



where a is the Goldstone field and 

oo _ d(t + a) d(t + a) nOl0 



dx a 



(2) 



Using the ADM formalism, we write the metric as ds 2 = 
-N 2 dt 2 + ^(dx l + N l dt)(dxi + N j dt), and expand at 
first order N — 1 + 4>, Ni — diB and 7^ = a 2 {(l - 
2ip)5ij + 2didjE}. Then the perturbative expansion of 



z 2 = 2 M 2 a 2 (e + n) 



(8) 



so that the quadratic action, after dropping more total 
derivative terms, reads 

= \J drd?x{{v'f cl{d iV f + ^V} , (9) 

where t is the conformal time and ' = ^ . The resulting 
mode equation for the Fourier transform v of v is 



v"+ (c 2 k 2 -^v = 0. 



(10) 



B. Modeling the non-adiabatic change of c s 

We are interested in the effects of a non-adiabatic 
change in c s originating from that of /1, i.e., of We 
will consider a small jump in fi during a short time about 
t = t . Since the overall normalization of the power 
spectrum is inversely proportional to c s , data allow only 
for a small jump |Ac s | < 1 in the speed of sound. For 
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simplicity we will assume that c s = 1 before the tran- 
sition. Since for c s close to unity (i.e., \i <C e) one has 
c s ~ 1 — we parametrize 



Mr) 



r , T < T , 

-2 Ac s • /(r) , r < t < r + At , (11) 
I -2 Ac s , t > T + At , 



with /(to) = and /(to + At) = 1. Here the magnitude 
of the jump, Ac s , has to be small in order to match 
observations, but for a sufficiently short At the rate of 
change can be large enough to open the possibility of an 
observable signature. For the sake of simplicity, one may 
be tempted to choose / to be linear in t, or even a step 
function. However, a discontinuity in c s or its derivatives 
can generate UV artifacts that contaminate the relevant 
information. For example, if fi jumps as a step function, 
then the feature coming from the sudden transition at 
t = to persists in all scales k > \to\~ 1 already at the 
two-point function level [15, 16]. With / linear in t, 
such that jU is continuous but pi is not, the oscillations in 
the power spectrum die off at large k, but the bispectrum 
is linearly divergent. It turns out that in order to obtain 
a vanishing bispectrum in the UV we must choose \x to 
be a function of class C 2 . 



C. Mode function 

In order to solve eq. (10) we employ the iterative 
method of [20] , to compute the mode function as a power 
series in the small parameter Ac s . We will denote by 
subscripts i, 2 and 3 the time intervals of before, dur- 
ing and after the period of non-adiabatic change in c s , 
respectively. 

Before the transition (r < To), eq. (10) becomes 



v'l + k 2 



2 + 3(2e-7?) 



) vi = 



0. 



(12) 



where r\ 



H 



is the second slow-roll 



2 H H ~ " 2«e 

parameter. Assuming e' to be negligible, so that rj ~ e, 
the solution of eq. (12) with the conventional boundary 



condition, v\(t — > —00) = e* fcr /v / 2fc, is 



Vl 



(13) 



Next, we introduce the dimcnsionless variables x = —kr 
(so that < x < +00, xq = —kro and Aa; = k At) and 
s = \f2kv. In terms of these variables eq. (13) reads 



si(x) 



(14) 



During the transition^ < t < t + At), we write the 
mode equation as 

d 2 s 2 / 2 + 3e 



dx 2 



(l -^^s 2 {x) = 2Ac s p{x)s 2 {x) + ... , (15) 



where ... denotes, from here on, any 0(e 2 , Ac 2 , e Ac s ) 
correction and 

Requiring that s 2 satisfies the boundary conditions 
^2(^0) = si(xo) and s 2 (xo) = s'^xo) fixes the homo- 
geneous solution of (15) to be just Si(x). Using the stan- 
dard Green's function technique, we obtain 



s 2 (x) = si(x) + i Ac s / dyp(y)s 1 (y) (17) 

J X 

x {s*i(y) si(x) - sl(x) Sl (y)} + ... . 
For t > tq + At, the mode equation reads 



d 2 s 3 



+ 



(l-2, 



2 + 3e 



da; 2 Vl-2Ac s a; 2 
whose solution is 



)« 3 (aO=0, (18) 



S3 



LhV . 

l i+ £ Vvl-2Ac s 



The coefficients a and (3 are determined by matching 
S3 and s 2 smoothly at x = xq — Ax: S3(a;o — Aa;) = 
s 2 (xq — Ax) and s' 3 (xq — Ax) = s' 2 (xq — Ax). 



D. Power spectrum 

To show concrete results, we choose our / to be the 
simplest polynomial in t that gives a C 2 interpolation 
between /(to) = and /(t + At) = 1: 



f(r) = 



(t - t ) 3 {3(2t - 2t - 3At)(t - t - At) + At 2 } 



At 5 



(20) 

Omitting the gruesome details of the intermediate steps, 
we eventually obtain 



a = 1 - - Ac s (2 a; - Aa;) + . . . 



(21) 



- t ^°- Ax) ^{^Ax-(l-^f)^} + 



2 Aa; 
giving the power spectrum 
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V(k) 



k " lim 



H 2 



2tt 2 x^q 2 k z(x) 2 
H 2 



8ir 2 M 2 e 



8ir 2 M 2 e 
{l + 2e (2- 7is -log2A;) 



{l - Ac s + 2e (2 - 7jE - log2fc) + ... } |c 



(22) 



-Ac, 



1 + 45 cos fc(2r + At) 



sin fcAr cos kAr sin kAr 



3fc 3 Ar 3 



fc 5 Ar 5 



}■ 



A/„- 



r P(k) 



1.35 






1.25 






1.20 






1.15 






-1 


1 2 


3 ^ 



for this system. We will use here the fact that a suffi- 
ciently rapid change in the speed of sound will mostly af- 
fect the modes that are deep inside the horizon at r ~ tq. 
For the analysis of those modes it is sufficient to study 
the action (1) in the decoupling limit. 

Let us first check the regime of validity of the quadratic 
Lagrangian in the decoupling limit. In this limit, the 
quadratic part of the Goldstone sector of (1) is 



Sl 2) = / drd 3 x 



M 2 H 2 e 



<j' 2 -c 2 (d i a) 2 -3U 2 ea' 2 



(24) 

Log 10 ( kT a ) w ^ n = a < j a Upon canonical normalization of a by 
w = -^f-cr, we get 



FIG. 1. Spectral density for e = .01, Ac s = -.005 and At = 
— .01 ro. Solid blue line: eq. (22), dotted red: IR asymptotics, 
dashed red: UV asymptotics. 



where is the Euler-Mascheroni constant. The asymp- 
totic behavior of the two-point function 



V(k) 



H 2 
T 2 M„ 2 e 



[1 + 2e (2 - lE 



[1 - Ac s + 2e (2 



log2fc)] , 
\kr \ «1, 

- 7s - log 2A:)] 
\kr \ > 1, 



(23) 

agrees, as expected, with the spectral densities for c s — 1 
and c s = l+Ac s +0(Ac 2 ) without any sudden transition. 

The feature generated by the non-adiabatic change in 
c s shows up in the term proportional to Ac s in (22): it 
vanishes for small fc, as can be seen from the IR asymp- 
totics. Oscillatory behavior of the amplitude of Ac s ap- 
pears once we deviate from the IR regime, i.e., —k tq > 1, 
but it dies off at k ~ 1/At due to the negative powers of 
k At in the oscillating terms. A plot of V(k) for e = .01, 
Ac s = —.005 and At = —.01 tq is shown in Fig.l. 



III. DECOUPLING LIMIT AND BISPECTRA 

In the previous section we have derived the two-point 
function of the perturbations by computing the full gauge 
invariant quadratic action with a varying speed of sound. 
In order to compute the three-point function of the per- 
turbations we need to move on to the cubic Lagrangian 



w' 2 - c 2 (d lW f + —w 2 



2He^- +eH 2 (3c 2 s + l-2r 1 ))w 2 



, (25) 



where the second line contains terms that should be ne- 
glected in the decoupling limit k ^> y/eH [17]. 

In particular, in the limit At <^ |t | the most im- 
portant difference between (25) and the exact action (9) 
comes from the term proportional to z 1 / z. Therefore the 
analysis of the action (24) will be consistent with the 
exact one of section II as long as z'/z ~ Ac s /At < H 
(where we have used c s ~ 1), i.e., 



I Ac J < 



At 

R 



(26) 



We expect that the regime of validity of the decoupling 
limit at the quadratic level will correspond to that at the 
cubic level. In other words, if the cubic Lagrangian in 
the decoupling limit turns out to be non vanishing, then 
it will give the dominant contribution to the three-point 
function up to 0(y/e) corrections. This procedure has 
been shown to give correct results, e.g., in the regime 
of small speed of sound [17] and has been invoked for 
the study of models with resonant nongaussianities [21]. 
While a formal proof of the validity of the decoupling 
approximation for the cubic Lagrangian is beyond the 
scope of the present paper, we believe that the explicit 
check we have provided for the quadratic Lagrangian, 
together with the examples described above, supports 
the use of this approach. 
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Let us then move on to the cubic action: 



be valid if sff < S y /> , with 



j(2) 



^ 3 > = Jdr d 3 x ^ { r/V acx' 2 + H 2 v a' 3 - H 2 n a'(d t a) 2 
+0(eH 2 )} = J dr(-H int ), (27) 

and the bispectrum is given by 

{n(k U T)n{k 2 ,T)n{k 3 ,T)) (28) 

= -iH 3 I dT'([a(k U T)a(k 2 ,T)a(h,r), H int (/)]>, 

where we have used 1Z = —Ha. 

Before we go ahead and compute (28) we must make 
sure that we trust our perturbative expansion. The most 
relevant term in the cubic action is the term proportional 

?(3) 



c(3) 



a Afi c 2 Ha 
{e + n) a ~ At" T ~H~ 



~ 10- 5 |Ac, 



Ar ' 
(29) 

where in the last equality we have used the fact that 
c s ~ 1, and where Ha ~ 10~ 5 is the amplitude of the 
perturbations. Let us also note that At/|to| is the dura- 
tion of the process in units of Hubble time. We thus find 
that as long as » 10~ 5 |Ac s | we trust our perturba- 
tive expansion. Since this condition is less stringent than 
(26), (26) is the only condition we require for the validity 
of our analysis. 



A. Equilateral case 

We are now in position to evaluate the three-point 
function. We first consider the equilateral configuration. 



to fi , let us denote it by S , . Perturbation theory will With cr^ 
I 



Hz 



we have 



5M 2 H 4 k 6 
6tt 4 -P(fc) 2 



[-'{ 









50 
81 



Ac 



-21 



.20 At 5 



+ /° -((-n 3 -f^K) 2 )}- 

Jto+At t V 2 >> 

1 ; /to At (3 to -4 At) 8t -7At\ 

P + fc3 J cos 3fcr ° 



0-3(1" -> 0) 3 + c.c. 



( At (t + At) (3 t + 7 At) 8to + 15At\ 
+ ( — p J cos3fc(T + At) 



-( 



/c 2 

3 t 2 At 2 12 t 2 - 60 t At + 1 1 At 2 22 

2k 6k 3 + 3fc5, 

3At 2 (t + At) 2 12 t 2 + 84 t At + 83 At 2 22 x v. 

+ 3~F J sm3fc ( T ° + Ar vJ ' 



2k 



6fc 3 



sin 3/cto 

22 > 
3fc5- 



(30) 



It can easily be found that the maximum of the feature 
is 

/N q L , max -l-7|A Cs |^{l + 0(^)}, (31) 
and occurs at 

*™~ = S{ 1+ °(f)}' < 32 > 

A plot of for Ac s = —.005 and At = —.01 t is shown 
in Fig. 2. The beats appearing in Fig. 2 originate from the 
fact that the function /x(r) is not infinitely differentiablc. 
We have checked numerically that the use of a smooth 
fi(r) cx (tanh ( 1 ^ L ) + l) leads to the disappearance of 
the beats without modifying the qualitative behavior of 



the three-point function. In particular, the scalings (31) 
and (32) found above are unaffected (up to an overall 
O(l) factor) by the choice of a C 2 -function /x(t). 

Note that the bispectrum (30) does not vanish in the 
limit k — > 00, that is, c s 7^ 1 leads to nonvanishing non- 
gaussianity even in the decoupling limit. The limiting 
value agrees with that found in [22, 23] for constant c s . 



B. Squeezed limit 

In the squeezed limit (fci = k 2 = k% = k$), the 

contribution to /nl from the jj! a a' 2 interaction is 
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75 Ac s f 2 / T Ar 2 
4 Ar 5 I V 3 fc? 



/At(2t 



At) 



2kl 



TQ 



2Ar 



sin 



COS 2fcLT — f 

At (2t 



At 2 (t + At) t + 3 At 



3k 2 L 



3At) 



2fc£ 



COs2fci(To 

sin2A: L (To + At)| 



At) 



O ^ 



(33) 



50 



-50 




Log I0 (-*T ) 



FIG. 2. Equilateral / NL for Ac s = -.005 and At = -.01 r . 



Its maximum, 



'NL.max 



= 1.3 I At 



l^ol 
At 



1 + 



At 

TO 



occurs at 



2J5 

a7 



1 + 



TO 



(34) 



(35) 



Note that it is smaller than that of the equilateral case - 
as expected, since the amplification of the a modes occurs 
deep into the horizon - by a factor of At/|to|. The lead- 
ing contributions to /nl from the other interactions are 
suppressed by ©(fcg/fci), and therefore negligible. Us- 
ing (33) and (22), one can also check that the consistency 
relation of [19] 



/nl(^£o ks 



5 dlnP(fc) 



12 dlnfc 
holds up to O^Acg, corrections. 



(36) 



General case 



Figure 2 shows that the behavior of f^{k) is given 
by high frequency oscillations modulated by an envelope 
function. We have checked that the envelope of can 
be written as 



.4 



eq 



B cq + C oq sin(3fcAT + /3 eq ) , 



(37) 



I 

with A eq , J5 eq , C cq and /3 eq functions of k, To and At. 
Since we factored out the high frequency parts, e 3kT ° , 
obtain (37), may be written as 



to 



J] 



NL 



A eq +sin(3/cTo- 



B r 



Cn, 



y eq sin(3fcAT + /3 eq ) . 

(38) 

we get the similar result with 3 k 



In the squeezed limit, 
replaced by 2k l- 

The above results can be generalized to arbitrary con- 
figurations of k\ , &2 and k% . The complete expression for 
fNh(ki,k 2 ,k 3 ) is too lengthy to be explicitly presented, 
and it can be schematically written as 



fNL(ki,k 2 ,k 3 ) = A+sm(KT +a)y/B + Csm(KAr + /3) , 

where K — k\ + ki + ^3, and A, B, C, a and /3 are slowly 
varying functions of k\, fc 2 , k 3 , t and At. This might be 
useful as a template for data analysis. We also checked 
that the At — > limit of (39) took the form of 

fNL{kx,k 2 ,k 3 ) -> A + Bsin(KT O + 0), (40) 

At->0 

which agrees with the Ansatz proposed by, e.g., [24]. 



IV. DISCUSSION 

We have shown that a sudden change in the speed of 
sound generates various distinct features in the primor- 
dial spectrum of scalar perturbations. 

First, due to the fact that the spectrum of perturba- 
tions is inversely proportional to the speed of sound, the 
normalization of the two-point function for modes that 
were deep inside the horizon before the transition took 
place is offset by an amount ~ 1 — Ac s with respect to 
that of the modes that were already far outside the hori- 
zon at the same instant. 

Second, the two-point function of the scalar perturba- 
tions shows an oscillating pattern of amplitude Ac s that 
dies off at large k. But oscillations last all the way to 
k — > oo for an infinitely sharp jump - in agreement with 
the results of [15, 16]. 

Third, the three-point function (in the equilateral con- 
figuration) shows a feature of maximum amplitude ~ 
Ac s (to/At) at scales k ~ At -1 . In particular, it is 
quadratically divergent at large k in the case of an in- 
stantaneous transition (At = 0). An analogous behav- 
ior has been recently found in [6], who have studied a 
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rapid transition in the inflationary potential and have 
shown that the seemingly violent UV divergence in mo- 
mentum space corresponds, due to its oscillatory nature, 
to a much milder logarithmic divergence in real space. 

The three-point correlation function found in the 
present paper shares many properties with that found 
in [6]. As a consequence, the discussion of the detection 
prospects for such a three-point function given in [6] ap- 
plies to our scenario as well: if the scale 1 / 1 To | roughly 
corresponds to current size of the horizon and the dura- 
tion At of the transition is short enough, a feature like 
(30) or (33) in the CMB bispcctrum might be detectable 
in the upcoming surveys. 

Let us conclude by emphasizing that the parameter 
controlling the speed of sound is just one of many pa- 



rameters appearing in the effective field theoretical de- 
scription of the inflationary perturbations [17]. It would 
be interesting to investigate whether a sudden change in 
some other parameters (e.g., those controlling the behav- 
ior of the tensor modes) would result in novel observable 
features in the CMB or in the primordial spectrum of 
gravitational waves. 
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